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ABSTRACT 

We investigate the influence of the geometry of the solar filament magnetic structure on the large- 
amplitude longitudinal oscillations. A representative filament flux tube is modeled as composed of 
a cool thread centered in a dipped part with hot coronal regions on either side. We have found 
the normal modes of the system, and establish that the observed longitudinal oscillations are well 
described with the fundamental mode. For small and intermediate curvature radii and moderate to 
large density contrast between the prominence and the corona, the main restoring force is the solar 
gravity. In this full wave description of the oscillation a simple expression for the oscillation frequencies 
is derived in which the pressure-driven term introduces a small correction. We have also found that 
the normal modes are almost independent of the geometry of the hot regions of the tube. We conclude 
that observed large-amplitude longitudinal oscillations are driven by the projected gravity along the 
flux tubes, and are strongly influenced by the curvature of the dips of the magnetic field in which the 
threads reside. 



1. INTRODUCTION 

Large-amplitude longitudinal (LA L) oscillat i ons i n 
prominences were first reported by IJing et alJ (|2003h : 
since then, only fe w additional reports of these mo- 
tions have appeare d () Jing et al.ll2006l ; IVrsnak et al.ll2007t 
IZhang et al.l 12012] ). These oscillations produce motions 
along the magnetic field, have long periods of 50-160 min- 
utes, and are damped in 2.3-6.2 cycles, with high velocity 
amplitudes in the range 30-100 kms -1 . LAL oscillations 
are apparently triggered by an energetic event: a sub- 
flare, a microflare, or a flare close to the filament. 

Several models have been proposed to explain the 
restoring force and da mping mechanism of the LAL oscil- 
lations (see review bv lTripathi et al.ll2009h . but most do 
not successfully describe the thread motions. Recently, 
we studied the oscillations of th reads forming the ba- 
sic components of a pro minence dLuna &: Karpenll2012f) 
in a 3D sheared arcade ()DeVore et al.ll2005t iLuna et al.1 
|2012| ). We found that the restoring force is mainly the 
gravity and the pressure forces are small. This type of 
oscillation resembles the motion of a gravity-driven pen- 
dulum, where the frequency only depends on the solar 
gravity and the flux-tube dip curvature. We estimated 
the minimum value of the magnetic field at the tube dips 
and found agreement with previous estimates and ob- 
served values. Additionally, this study revealed a new 
method for measuring the radiu s of curvature of the fil- 
ament dips. IZhang et al.1 (|2012f) observed and analyzed 
an oscillating prominence and found that the motion is 
prod uced along a dipped m agnetic field, in agreement 
with ILuna fc Karpenl il2l)12i. These studies reveal that 
the LAL oscillations are strongly related to the filament- 
channel geometry. 
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The filament- channel structure is not well understood 
and several models have been suggested. The sheared 
arcade and the flux rope models are the most successful 
candidates explaini ng most of the obser vational evidence 
(see the review by Mackay et a.1.1 [20Toh . In these mod- 
els the magnetic structure is static and independent of 
the prominence evolution because the plasma-/? is small, 
and the structure has dips where the cool prominence 
resides. Due to the low plasma however, most mod- 
els agree that these dips are not caused by weight of the 
prominence. The prominence mass forms in the dipped 
part of the magnetic field because it is a gravity poten- 
tial well where evaporated mass plasma conden se and col- 
lect (jAntiochos et al.lll999HKarpen et al.ll2003[) . There is 
also direct observational evidence from polarimctric in- 
yersions of the dippe d ma gnetic stru cture of the filaments 
(|L6pez Ariste et al.ll2006l : IXu et al.ll2012f ). In these mod- 
els the curvature in the dips is large and the LAL os- 
cillations could be strongly influenced by these geome- 
tries. In contrast, where the magnetic field is slightly 
curved by the prominence weig ht, the curvature effects 
on the oscillations are negligible (|Qliver et alJ|1992 l ll993l : 
lOliver fc Ballesterlll995l : iTerradas et al.H2001l L ~ 

In our previous study of LAL oscillations, we as- 
sumed the threads to be solid masses moving in curved 
flux-tubes without interaction with the surrounding hot 
plasma. In this work we use a full wave description of the 
oscillation, expanding our previous investigations. We 
focus on the restoring forces of the LAL oscillations and 
the influence of the curvature of the filament magnetic 
fields in different tube geometries, and compare the re- 
sulting thread motions with observed LAL oscillations 
properties. 

2. FLUX TUBES WITH CURVATURE 

In this work we assume that the plasma is low-/3 and 
confined with static magnetic field. In this regime the 
plasma motion is d escribed with Equations (l)-(4) of 
iKarpen et al.1 ()2005l ). We additionally consider that the 
system is adiabatic, with no heating and radiation, and 



2 



M. Luna et al. 



the tubes have a uniform width. Thus, the terms associ- 
ated with the energy loss and gains, and the area expan- 
sion, can be neglected. We linearize this set of equations 
to obtain the equations for the perturbed quantities 
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where the "0" index means the equilibrium quantity that 
depends on the coordinate s, the "1" index means the 
perturbed quantity, v is the perturbed velocity (a zero 
background velocity is considered), and gu is the gravity 
projected along the flux tube. Additionally the plasma 
is in hydrostatic equilibrium 
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Combining Equations ([I])-® we obtain the equation for 
the velocity perturbation 
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where c s — c s (s) — 7Pn( s) / Pais). This equation is s imi- 
lar to Equation (7.30) of IGoedbloed fc Poedtsl (|200l as- 
suming a displacement constrained along the flux tube. 

We model the filament flux-tube geometry as com- 
posed of up to 3 curved segments (see Fig. [T]). These 
segments are contained in a vertical plane with the cen- 
ter of curvature located above or below the tube for a 
concave-up or concave-down segments respectively. Each 
segment has a constant radius of curvature R. The ra- 
dius of curvature is positive for a concave-up segment 
and negative for a concave-down segment. With these 
considerations, on each segment of the tube Equation 
© is 
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where l/R = d6/ds defining 9 = s/R. Assuming now 
that the radius of curvature R is sufficiently large to fulfill 
the condition \s/R\ <C 1, the 9 angle takes small values 
in all the positions of the tube. Thus considering linear 
approximations of the sinusoidal functions (sin 9 = 9 and 
cos 9=1), then 
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This e quation reduces to Equation (8) of lDiaz fc Roberts! 
(2006) for a low-/? plasma when curvature is negligible 
(i.e., R — >• oo). Hence, for a straight tube the equation 
describing the perturbations is simply 
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We consider solutions for which the time dependence 
is a simple harmonic oscillation with frequency u> and 
the velocity perturbation takes the form v{s)e~' tult . Thus 



Equation ([7]) becomes 
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This equation describes the oscillatory motion of a 
plasma embedded in a curved tube with constant cur- 
vature in hydrostatic equilibrium. 

2.1. Isothermal plasma flux-tube segments 

Solving Equation © is complicated because the sound 
velocity depends on the position along the tube, i.e. c s = 
c s (s). However, if we assume that the segments of the 
tube are isothermal, the sound velocity is uniform on 
each segment, i.e. c s ^ c s (s). Thus, we can perform a 
linear change of variable defined by 
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and Equation [3] becomes 
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We consider that the plasma within the flux tubes of 
the filament is distributed in three isothermal regions. 
The thread consists of cool plasma centered at the dip 
with an uniform sound speed, c sp . The remaining ther- 
mal regions are the hot coronal plasma filling the tube 
from both ends of the thread to the footpoints, with a 
uniform sound speed, c sc . In this work we ignore the 
small transition region from the cool prominence thread 
to the coronal hot plasma (the so-called PCTR) at both 
ends of the thread. This region is thin in comparison 
with the thread length, and the influence on the thread 
oscillations may be small, so we leave the incorporation 
of this region for a future study. We note that both the 
equilibrium density and pressure depend on s, despite 
having a uniform sound speed along each flux-tube seg- 
ment. In this situation the density and pressure decrease 
exponentially with height, at a rate given by the pres- 
sure scale height. For a typical hot coronal plasma the 
pressure scale height is a few tens of Mm. However, for 
a typical prominence it is about 0.2 Mm. Therefore the 
pressure and density are roughly uniform in the coro- 
nal parts of our low-lying flux tubes, whereas they vary 
rapidly along the curved thread. 

Equation (jTTJ) governs the plasma motion on each part 
of the piecewise flux tube with uniform temperature and 
R. This equation is a second-order ordinary equation 
with a irregular singular point at infinity, so its solu- 
tion can be expressed as a combination of two regular 
linearly independent functions at each point (except at 
infinity). The general solution can be written in terms 
of c onfluent hypergeometric func tions (CHFs) M(a,b;x) 
(see I Abramowitz fc Stegunlll972D as 
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with A\ and A2 being two arbitrary constants and 

, 2 [Rlo 2 \ , 

A = ~ 1 . 13 

7 V 9o J 

This solution can be expressed in terms of the well-known 
Hermite polynomials if A = 2n, with n being an integer, 
but even in this case there is a second regular solution 
which cannot be discarded. Hence, it is suitable to work 
directly with the general solution in terms of CHFs (Eq. 

Ull). 

Regarding the symmetry at r = 0, the CHFs have the 
following property at the origin: 



M(a,b,0) = 1, 



(14) 



so we obtain a symmetric function by setting A2 = and 
an antisymmetric one by setting A\ = 0. 

In this work we use three flux-tube models to inves- 
tigate the influence of the hot region on the motion of 
the cool plasma. The cool thread is clearly located in a 
concave-u p field depression or dip as shown in our ear- 
lier study (|Luna et al.ll2012l ). Model 1 is the simplest: a 
tube with straight hot regions connecting to the chromo- 
sphere (see upper sketch of Fig. [TJ. Model 2 is a tube 
which is concave-up in all domains (see middle sketch 
of Fig. []J>). Finally, Model 3 has the most complex ge- 
ometry considered in this work, as shown in the bottom 
sketch of Figure [TJ a dipped part containing a central 
cool thread and two hot plasma regions at both sides, 
connected with two concave-down segments forming the 
legs of a M-shaped tube. 




Fig. 1. — Sketch of the configurations considered in this study, 
differentiated by the shape of the flux tube in the hot plasma region. 
The upper sketch corresponds to Model 1, the middle sketch to 
Model 2, and the lower sketch to the M-shaped Model 3. The 
length of the flux tubes is 2L, the length of thread is 21, and the 
dipped part of the tube has a length 2d in Model 3. 



2.2. Boundary conditions 

The solutions of each region of our piecewise flux tube 
model must be joined by the boundary conditions. The 
resulting solutions are the normal modes. To solve this 
problem we need two types of boundary conditions: line- 
tying at the chromospheric footpoints of the tube (lo- 
cated at s = ±L, with 2L being the length of the sup- 
porting magnetic flux tube) and jump conditions at the 
plasma interfaces. We adopt the simplest possible con- 
ditions at the fo otpoints: a line-tied rigid w all with no 
flow (|Hoodlll986t Ivan der Linden et al.lfl99l IDi'az et al.l 



120041) . namely 



v(±L) = . 



(15) 



appropriate mainly because there is not enough energy in 
the perturbations carried by the coronal material to drive 
important motions in the chromosphere. Moreover, this 
simple condition completely decouples the coronal and 
the photospheric plasmas (so the dynamics of the photo- 
spheric plasma are not taken into account), although we 
plan to address more advanced "flow-through" conditions 
and coupling with the chromosphere in future works. 

The second set of boundary conditions is prescribed at 
the interface between the cool and hot plasmas, located 
at s = I, with 21 being the length of the dense thread. We 
need to carefully deduc e these boundary conditions. Fol- 
lowing [Chandrasckhar (1961J we integrate the differen- 
tial Equations (p}- (J3j) across the boundary (from s = l — a 
to s = I + a) and let a — ¥ 0. The integral of any vari- 
able that does not have an infinite jump becomes zero 
in that limit. First, integrating the hydrostatic equilib- 
rium (Eq. [4j> we obtain [po] = with the commonly used 
notation [a] =0,2 — 0,1 for the jump between mediums 1 
and 2. Therefore the equilibrium pressure is continuous 
at the boundary. Similarly, from Equations CO)-© the 
two boundary conditions for the velocity are 



0. 



dv 
ds 



= 
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Since the chromosphere is much denser than the corona, 
it behaves as a purely reflecting layer. This condition is 



Thus, including gravity does not affect the jump 
conditions, in agreem ent with iDiaz et all (|2004l ) and 
IDi'az fe Roberts! pool) . 

3. UNCOUPLED THREAD 

In a first approximation we consider the oscillation of 
the thread alone, assuming that the hot plasma filling 
the rest of the tube is irrelevant to the dynamics of the 
prominen ce. This assumption i s justified by our previ- 
ous work (jLuna fc K arpen 2012), in which we found that 
the dynamics of the threads are basically governed by the 
gravity, and the interaction with the ambient hot plasma 
is small. In this situation the boundary conditions of 
Equation (fl6| are not applicable. We assume that the 
thread moves freely and thus the boundaries at s = ±1 
are open. The motion of the thread is described by Equa- 
tion (fl"2"|) ; the Ai, A2, and A parameters can be chosen 
freely because there is no constraining condition. 

In Figure [2j we have plotted different symmetric 
and antisymmetric solutions of Equation (fl"2|) along the 
thread, \s\ < I, for different A values. For a symmetric so- 
lution and A = the whole thread moves as a rigid body 
with an uniform velocity. For a small value (A = 0.3) the 
symmetric velocity distribution differs only slightly from 
the uniform case: the velocity at the ends of the thread 
is slightly smaller than in the center. For A = 10 the ve- 
locity clearly is not uniform but still positive. For A = 40 
the motion is more complex, however to with the thread 
ends moving out of phase with respect the center of the 
thread. Hence, the A parameter gives the spatial coher- 
ence of the velocity in the thread. For small A the spatial 
coherence is high, with a uniform velocity distribution, 
whereas for larger values the spatial coherence is small 
with steep velocity gradients along the thread. For anti- 
symmetric solutions and A = both ends of the threads 
move in opposite directions, compressing and rarefying 
the thread plasma, whereas the center remains at rest. 
For A = 0.3 the antisymmetric solution is slightly differ- 
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Fig. 2. — The symmetric (solid l ines ) and antisymmetric (dashed 
lines) normalized solutions of Eq. fl I L'li for the values A = (black), 
A = 0.3 (red), A = 10 (green), and A = 40 (blue) for a thread of 



half-length I = 5 Mm, R = 75 Mm, and &>, 



20 kms~ 



10 



ent, whereas the solution is more complicated for A 
and 40. 

The general solution is a superposition of these modes, 
where the exact combination depends on the initial con- 
ditions. One special case is worth studying carefully: an 
exact rigid-body perturbation in which the whole thread 
is shifted from the initial equilibrium position. In this 
case the solution is symmetric and the spatial coherence 
parameter is A = 0. Thus the frequency of the oscillation 



This solution was described in iLuna &: Karpenl ((20121 
Eq. 4), and predicts the thread oscillating as a gravity- 
driven pendulum. For a general initial perturbation 
Ao(s,t = 0) this solution is no longer valid and Equa- 
tion (0 must be solved instead. 

4. CURVED THREAD WITH A STRAIGHT FIELD IN THE 
HOT PLASMA REGION 

We consider now Model 1, where the thread is located 
in the center of the flux-tube dip and the remainder of the 
flux tube is straight and filled with hot coronal plasma 
(upper sketch of Fig. [1} . The solution to this piecewise 
model is given by Equation (|12[) in the dense region and 
Equation (JS) with go/R — > in the evacuated region, 
namely 

v(s) = { + B 2 r P M (-^, §;r-2) , < a < I, (18) 
D x sin j^(Z, - s)\ , l<s<L. 

We have ensured that the line-tying boundary condition 
(Eq. II 5|) is enforced by choosing only the sine function in 
the hot region. Using the series expansion for the CHF, 
it can be proved that 



/R^o V 2 19o/R ' 2 ' 2c s 2 R 
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so we recove r the fully-s t raigh t ened case with very large 
R studied in iDi'az et all ((20021 IMol ). 



Next, we apply the boundary conditions (Eq. [TB]) at 
the interface between the cool and hot plasma at s = 
±1. We obtain the following dispersion relation for the 
symmetric normal modes (symmetric with respect to s = 
0) after eliminating the amplitude constants B\ and D\, 
and considering B2 = 0: 
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We define a new set of dimensionless variables as 
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With these definitions the dispersion relation takes the 
form 
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The parameter x is usually called density contrast when 
the density and pressure are assumed uniform in each 
segment along the tube. However, in this work we con- 
sider nonuniform densities and pressures. The variable 
X is equal to the density contrast only at the interface 
s = I because the pressure must be continuous across the 
thread-corona interface (see i J2.2l) . Thus x — Pp(l)/Pc(l) 
and we call this parameter the contrast or density con- 
trast hereafter. 
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Fig. 3. — The symmetric (solid lines) and antisymmetric (dashed 
lines) normal modes for Model 1 as a function of the radius of 
curvature R of the dip. The fundamental mode is plotted with 
a t hick solid line. The gravity-driven frequency, u s , defined by 
Eq. 1171 (dot-dashed line), and the pressure-driven frequency, w s , of 
Eq. 1261 (three dot dashed line), are also plott ed. Additionally, the 
approximate fundamental frequency of Eq. 1271 is shown (dotted 
line). The half-length of the tube is set to L = 100 Mm, the thread 
half-length I = 5 Mm, the coronal sound speed c sc = 200 kms -1 , 
and the contrast \ = 100. 

Similarly, we obtain the dispersion relation for the anti- 
symmetric modes by eliminating the amplitude constants 
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B2 and Di, and considering B\ = 
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in terms of the dimensionless variables of Equation ([2~Tj) . 

In Figure [3] we have plotted the fundamental mode 
and several overtones of the dispersion relations (Eqs. 
|2"21 and [53]) as a function of the dip radius of curvature 
R. The fundamental mode (solid line) decreases with 
R, clearly indicating that the fundamental mode is af- 
fected by the dip curvature. In contrast, the frequencies 
of the overtones are independent of R, indicating that 
these overtones are purely sound-like or pressure-driven 
modes. We have also plotted the curve u> = uj g (or equiv- 
alently — 1), which corresponds to the case where the 
restoring force of the thread oscil lation is exclusively th e 
gravity projected along the tube ([Luna &: Karpenl l2012). 
We see that the fundamental mode is very similar to the 
gravity-driven oscillation oj w uj g for small values of R. 
For an intermediate radius of curvature, R — 200 Mm, 
the gravity-driven approximation differs by 18% with re- 
spect to the exact solution, and by 38% for a relatively 
large radius R = 600 Mm. Therefore the approximation 
uj ~ cj g is good for small and intermediate values of R. 
The increasing deviation from the purely gravity-driven 
case indicates that the pressure force has some influence 
for intermediate values of R, and becomes important for 
large radii. 

The argument in the CHFs of the symmetric dispersion 
relation (Eq. [22]) is small (x-fW 2 ® 2 /2 < 1) for the values 
of the parameters considered here. Additionally, 57$ = 
luL/c sc < 1 for the fundamental mode. Thus, we expand 
the dispersion relation and obtain 
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Substituting the dimensional variables (Eq. |2T]) yields 
an approximate expression for the fundamental mode 
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The fundamental frequency has two contributions: the 
gravity-driven frequency uj 9 , and the pressure-driven 
slow oscillation frequency 



l(L-l) X 
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Now Equation ([25]) can be written as 

2 22 

w fund = w 9 + w s • 

Figure (J3J exhibits a perfect match between the ap- 
proximate expression and the exact values of the funda- 
mental mode frequencies for a range of R values. The 
w g and uj s terms are identical at R w 330 Mm, while for 
larger values of R the major contribution to the restor- 
ing forces is the pressure gradient. However, the funda- 
mental mode differs significantly from uj s , except at huge 



values of R. lOliver et all (|1993ft and lOliver fc Ballesterl 
(1995) studied the oscillation modes of a slab in a 
Kippenhahn-Schliiter magnetic configuration with a very 
slight curvature. They concluded that the forces respon- 
sible for the slow modes are the pressure gradients, and 
the frequency of the fundamental mode is well described 
by uj s , consistent with the present result. 

Similarly to the symmetric case we find an approximate 
expression for the frequency of the first overtone. In the 
range of small I and large R the term W§ is small in the 
antisymmetric dispersion relation (Eq. I23p . The only 
way to balance the first term of the right-hand side of the 
equation is for the argument of the cotangent function to 
be near the first asymptote; then 



(1 - w)$ ' 

In dimensional variables this equation is simply 



L-r 



(28) 



(29) 



which is the frequency of a standing wave trapped in one 
of the hot sections of the tube with wavelength 2(L — I). 
In Figure [3] we see that the frequency is more or less 
independent of R. 

In Figure [4] the velocity perturbations of the first three 
normal modes are shown. In the fundamental mode the 
maximum velocity is centered within the thread and the 
thread speed is quite uniform, indicating that motion 
resembles the motion of a solid body. In the coronal 
part the velocity reduces to zero at the footpoints as 
dictated by the boundary conditions. The first overtone 
is an antisymmetric mode with a node (v = 0) at the 
center of the tube. In this mode the motion consists of 
compressions and rarefactions of the cool plasma with 
no net displacement of the thread. The second overtone 
has a node at each thread end, so the motion of the 
plasma is complex in this mode. All overtones reach their 
maximum velocities in the coronal parts of the tube. 
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Fig. 4. — The normalized velocity along the flux tube, v (s), for 
the fundamental mode (solid line), the first antisymmetric overtone 
(dashed line), and the first symmetric overtone (dot-dashed lines). 
The flux tube half-length is L = 100 Mm, the thread half-length 
is I = 5 Mm, dip radius of R = 75 Mm, coronal sound speed 
c sc = 200 kms" 1 , and contrast \ = 100. 

The frequencies of the symmetric and antisymmetric 
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normal modes as a function of the contrast \ are shown 
in Figure [5] The fundamental mode rapidly reaches a 
constant value for increasing values of Xj reflecting a 
weak dependence on the contrast. The frequency of the 
approximation Wf un( j is very similar to the fundamental 
mode. For relatively large values of the contrast the 
match between both curves is very good. The frequency 
lu s is inversely proportional to x an d thus rapidly goes 
to zero (triple dot-dash line in Fig. [5|, so the main con- 
tribution to the fundamental mode (Eq. I27[) comes from 
the gravity term u g . Then, for sufficiently large values 
of x, the oscillation of the thread decouples from the hot 
plasma and is governed by the gravity forces as shown 
in fj3j The frequencies of the overtones have also a com- 
plex dependence on x, with avoided crossings between 
modes. These avo i ded c rossin gs are similar t o thos e 
found in IDiaz et all (|200l and iDi'az fc Roberts! (|2006l) . 
and have the same origin: two modes cannot have the 
same frequency, since it would violate the conditions in 
a Stiirm-Liouville problem such as our Equation (|TT|) . 
Furthermore, the modes can have extrema in the promi- 
nence region and the hot region, and an avoided crossing 
occurs when the amplitude of the extrema in the promi- 
nence region changes from being small (such as the first 
overtone in Fig. 2]) to being larger than the amplitude in 
the hot region. 



range of large curvature radii. As in the previous section 
(Sj4]) we have found the dispersion relation for the normal 
modes of this model. Figure [6] shows that the fundamen- 
tal frequency is essentially identical to that of Model 1. 
The frequencies and spatial distribution of the velocity 
obtained with both symmetric and antisymmetric dis- 
persion relations are very similar to the corresponding 
solutions of Model 1 (0 for all values of R, L, x, and I. 
For this reason we do not show the dispersion relations 
for Model 2. 
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Fig. 5. — Same as Figure \3\ for the symmetric (solid lines) and 
antisymmetric (dashed lines) normal modes for Model 1 as a func- 
tion of the density contrast \ at the interface between the cool and 
hot plasma. The fundamental mode is plotted with a thick solid 
line. The half-length of the tube is set to L = 100 Mm, the thread 
half-length 1 = 5 Mm, a coronal sound speed c sc = 200 kms -1 , 
and the radius of curvature R = 75 Mm. 



5. INFLUENCE OF CURVATURE IN THE HOT REGIONS 

We consider next the effects of curvature in the hot 
parts of the tube, which were assumed to be straight in 
the previous section. Two different cases are discussed. 

5.1. Model 2: concave-up field in the hot region 

We consider a tube with constant curvature in the 
entire domain (see middle sketch of Fig. [TJ. We as- 
sume that isothermal cool plasma occupies the region 
\s\ < I, while the rest of the tube is filled with isothermal 
hot plasma. In this model we have the limitation that 
\R\ 3> L in order to satisfy the small-angle approxima- 
tion \s/R\ <C 1 (see fj2|). Thus we can only consider the 



Fig. 6. — Same as Figure \3\ for the frequency of the fundamental 
mode as a function of R, for Model 1 (dotted line), Model 2 (solid 
line), and Model 3 (dashed line). In Model 3 we have condisered 
a dip half-length d = 20 Mm, and a radius of curvature of the legs 
R 2 = -100 Mm. 



5.2. Model 3: M-shaped flux tube 

Finally, we consider a more realistic model composed of 
three different regions (see bottom sketch of Fig. [T]). The 
dipped part is equivalent to Model 2, with the cool and 
dense thread of length 2d centered and filling a portion 
21, and the remaining portion of the dip is filled with 
hot rarefied plasma. This dipped part is connected to 
the chromosphere by two concave-down segments filled 
with hot plasma. The resulting flux tube has an "M" 
sh ape, similar t o the thread-bearing flux tubes identified 
bv lLuna et all (|2012] ). 

The new feature of this M-shaped tube is the concave- 
down region, where the radius of curvature is negative. 
Thus in this region the linear substitution in Equation 
(fTU|) must be replaced by 



\2c 2 R 2 



(30) 



with i?2 < 0, the radius of curvature of the concave-down 
region. With this new definition Equation ([9]) becomes 
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The solution of this equation can be written again in 
terms of CHFs as 
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This functional form must be used for d < \s\ < L, while 
for \s\ < d, Equation (fT2|) is still valid. Using the bound- 
ary conditions to match the solutions at s = I, s = d 
(Eq. I16|) and the line-tying boundary condition at s = L 
(Eq. I15[) . we obtain a very cumbersome dispersion re- 
lation which involves products of five CHFs. For this 
reason, this dispersion relation is not shown in here. 

Figure |6] shows that the fundamental frequency for 
Model 3 is very similar to that of Model 1. In this 
figure we have considered a radius of curvature in the 
flux-tube legs of |i?2| = 100 Mm (note that this is the 
smallest possible value allowed by the small angle con- 
dition, \R 2 \ > L). For larger values of |i? 2 | the differ- 
ence between the frequencies of Models 1 and 3 are even 
smaller, and zero when \R 2 \ oo. We have also studied 
the dependence of the frequencies on the other parame- 
ters of the system and found identical results. 

6. THE RESTORING FORCE 

In our theoretical model of prominence oscillation there 
are two restoring forces: the gas pressure and the gravity. 
Their relative importance is determined by the ratio 



R 



(33) 



where i?i; m = l(L — 0x5o/ c sc i s a reference radius of 
curvature determined by measurable filament properties. 
The restoring force of the LAL oscillation is mainly the 
solar gravity when this ratio is small, i.e., R <C i?ii m - 
In contrast, the oscillation is pressure-driven for a large 
value of the ratio or R Ru m . With the data shown in 
Figure[3]we find that Ru m — 325 Mm, so for smaller radii 
of curvature the gravity is the restoring force. The length 
of the threads and density contrast have a wide range of 
measured or estimated v alues (see, e.g., lLabrosse et al.l 
120101 : iMackay et al.l l2010h . Considering larger values of 
the thread length and contrast of Z = 10 Mm and x = 
200, we find Ru m = 1200 Mm, indicating t hat the grav- 
ity do minates even for almost fiat tubes. iZhang et al.l 
(2012) studied the observed longitudinal oscillations of a 
prominence, and found a small radius of curvature less 
than 100 Mm. Therefore the restoring force is certainly 
the gravity, as we would expect in those structures. Ad- 
ditionally, the relative importance of the two restoring 
forces also depends on the geometrical parameter R/L 
(Eq. [33]), which can be determined from the global ge- 
ometry of the filament magnetic structure. For the deep- 
dip fl ux tubes of a double sheared ar cade dDeVore et all 
I2005D . this factor is always R/L< \ (|Luna et al.ll2012l ). 

In |21 we discussed that the spatial coherence of the 
motion of the thread is determined by A. Combining 
Equation fp~3|) with Equation (f27|) , we find that in Models 
1 to 3 

A = -4- (34) 

Thus the relative importance of the pressure and grav- 
ity forces determines the way the thread moves. For 
a gravity-driven oscillation A = 0, the thread basically 
moves almost as a solid body. However, for a pressure- 



driven motion the thread compresses and rarefies, chang- 
ing its shape in the oscillation. For the reference values of 
Figure [3] and a radius of 75 Mm the coherence parameter 
is A = 0.3. The velocity profile for this parameter value 
plotted in Figure^ is quite flat in the thread, indicating 
that the motion is very similar to the displacement of a 
rigid solid body. 

7. DISCUSSION AND CONCLUSIONS 

In this work we have studied the influence of the cur- 
vature on the longitudinal oscillations of a prominence. 
We have considered three different models, in which the 
region where the thread resides is modeled as a tube seg- 
ment with uniform curvature R. The differences between 
the three models are in the hot regions of the tubes. 
We have found that the frequency of the fundamental 
mode is dependent on R and can be approximated by 
w fund = Wg+w s 2 , where w g and uj s are the frequencies of 
the gravity-driven pendulum and the fundamental slow 
mode of a straight tube respectively. For small and inter- 
mediate R the frequency is very close to oj g . We modeled 
a prominence with realistic dimensions and found most of 
the filament flux-tube s have dips with a radius of curva- 
ture around 75 Mm (jLuna et alj I2012t ILuna k, Karpenl 
120121 ). We also inferred the rad ius of curvature from 
the oscillation perio ds reported bv lJing et al.1 ()2003l ) and 
iVrsnak et all (I2007D. and found 152 Mm and 62 Mm 
respectively. IZhang et al.l (|2012D observationally deter- 
mined the radius of filament-dip curvatures to be less 
than 100 Mm. Thus, the observations and theoretical 
models are consistent with this range of small and inter- 
mediate R. For larger radii the pressure force becomes 
more important and the fundamental frequency differs 
from w g . The frequencies of the overtones are basically 
independent of the curvature of the tube, consistent with 
the slow nature of these modes. 

The fundamental mode is also weakly dependent on the 
density contrast, Xi f° r small contrast, but as x increases 
the fundamental mode rapidly reaches a constant value 
that coincides with Wf un( j or uj g (note that Wf un d ~ 
because u; s is very small for relatively large contrast). 
Therefore the oscillation of the thread decouples from the 
environment for relatively larger values of the contrast, 
and the thread oscillates with the frequency w g . Obser- 
vational estimates of prominence densities give a broad 
range of possible x values, but typ ically the density con - 
trast is 100 or larger (see review bv lLabrosse et aLl l2010'). 
Therefore the threads are in the range of large density 
contrasts. The overtones depend on x, indicating the 
sound-like nature of these modes. The spatial velocity 
distribution of the fundamental mode along the tube is 
symmetric with respect to the tube center; although the 
maximum is located at the thread center, the velocity in 
the thread is more or less uniform. In the fundamental 
mode the motion is mainly concentrated in the thread. 
The overtones produce complex compression and rarefac- 
tion motions with small or zero net displacements of the 
thread. 

In Models 2 and 3, with curved hot regions, the fre- 
quencies and the spatial distribution are very similar to 
the corresponding values for Model 1. We conclude, 
therefore, that the shape of the hot regions is irrele- 
vant to the longitudinal oscillation, and that using a 
straight field-line approximation in these regions gives 
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results that are accurate enough and much easier to com- 
pute. Hence, the approximation given in Equation (|27[) 
is quite robust, despite only being truly valid for Model 
1. The sound speed in the corona is very high, and the 
contribution of the first term in Equation ([7]) is larger 
than terms involving curvature. Thus, the resulting mo- 
tion in the coronal parts of the flux tube is well described 
by Equation ([8]) . 

The relative importance of the pressure and gravity 
forces is determined by the geometry of the dipped part 
of the filament flux tubes. We found that the oscillation 
is gravity driven when R <C i?i; m , where i?i; m is deter- 
mined by the properties of the filament flux tubes and 
the cool thread. For typical prominences Ru m is much 
larger than the radius of curvature of filament models or 
observationally inferred values. 

We conclude that the longitudinal oscillations of a 
prominence are strongly influenced by the curvature of 
the dipped magnetic fields. In the fundamental mode 
the gravity dominates for small and intermediate radius 
of curvature, consist ent with values in our multi-threaded 
prominence model (jLuna et al.l I2012T ) and with the val- 
ues derived from observed oscillations. Similarly, for rel- 
atively large density contrast between the prominence 
and the corona the main restoring force is also the grav- 



ity. Thus, the frequency of the LAL oscillations is given 
by w g = y^go/R, and the pressure forces introduce only 
a small correction, showing that the LAL oscillations are 
not pure slow modes. This demonstrates that the oj g 
expression is robust and can be used for seismology of 
prominences, as we showed in lLuna fc Karpenl J2012). 

In this work we have not s tudied the damp i ng of the 
observed LAL oscillations. In lLuna &: Karpenl (|2012f) we 
found that the damping is associated with mass accretion 
onto the threads and nonadiabatic effects. The former 
produces strong damping at the beginning of the oscil- 
lation, and the latter yields a weak damping through- 
out the oscillation. In order to have a full and self- 
consistent model of prominence oscillations we must per- 
form a nonlinear study, including the temporal variation 
of the prominence mass and the nonadiabatic effects. 
This will be the subject of a future work. 



This work has been supported by the NASA Helio- 
physics SR&T program. M.L. also acknowledges support 
from the University of Maryland at College Park and 
the people of CRESST. A.J.D. acknowledges the finan- 
cial support by the Spanish Ministry of Science through 
project AYA2010-18029. 



REFERENCES 



Abramowitz, M., & Stegun, I. A. 1972, Handbook of 

Mathematical Functions, ed. Abramowitz, M. & Stegun, I. A. 
Antiochos, S. K., MacNeice, P. J., Spicer, D. S., & Klimchuk, 

J. A. 1999, ApJ, 512, 985 
Chandrasekhar, S. 1961, Hydrodynamic and Hydromagnctic 

Stability (CUP, New York: Dover Publications Inc., 1932) 
DeVore, C. R., Antiochos, S. K., & Aulanier, G. 2005, ApJ, 629, 

1122 

Diaz, A. J., Oliver, R., & Ballester, J. L. 2002, ApJ, 580, 550 
Diaz, A. J., Oliver, R., & Ballester, J. L. 2010, ApJ, 725, 1742 
Diaz, A. J., Oliver, R., Ballester, J. L., & Roberts, B. 2004, 

Astronomy and Astrophysics, 424, 1055 
Diaz, A. J., & Roberts, B. 2006, Astronomy and Astrophysics, 

458, 975 

Goedbloed, J. P. H., & Poedts, S. 2004, Principles of 

Magnetohydrodynamics, ed. Goedbloed, J. P. H. & Poedts, S. 
(Cambridge University Press) 
Hood, A. W. 1986, Solar Physics (ISSN 0038-0938), 105, 307 
Jing, J., Lee, J., Spirock, T. J., & Wang, H. 2006, Sol. Phys., 236, 
97 

Jing, J., Lee, J., Spirock, T. J., et al. 2003, ApJ, 584, L103 
Karpen, J. T., Antiochos, S. K., Klimchuk, J. A., & MacNeice, 
P. J. 2003, ApJ, 593, 1187 



Karpen, J. T., Tanner, S. E. M., Antiochos, S. K., & DeVore, 

C. R. 2005, ApJ, 635, 1319 
Labrosse, N., Heinzel, P., Vial, J., et al. 2010, Space Sci. Rev., 

151, 243 

Lopez Ariste, A., Aulanier, G., Schmieder, B., & Sainz Dalda, A. 

2006, A&A, 456, 725 
Luna, M., & Karpen, J. 2012, The Astrophysical Journal, 750, LI 
Luna, M., Karpen, J. T., & Devore, C. R. 2012, The 

Astrophysical Journal, 746, 30 
Mackay, D. H., Karpen, J. T., Ballester, J. L., Schmieder, B., & 

Aulanier, G. 2010, Space Sci. Rev., 151, 333 
Oliver, R., & Ballester, J. L. 1995, Astrophysical Journal v. 448, 

448, 444 

Oliver, R., Ballester, J. L., Hood, A. W., & Priest, E. R. 1992, 
ApJ, 400, 369 

Oliver, R., Ballester, J. L., Hood, A. W., & Priest, E. R. 1993, 
ApJ, 409, 809 

Terradas, J., Oliver, R., & Ballester, J. L. 2001, A&A, 378, 635 
Tripathi, D., Isobe, H., & Jain, R. 2009, Space Sci. Rev., 149, 283 
van der Linden, R. A. M., Hood, A. W., & Goedbloed, J. P. 1994, 

Sol. Phys., 154, 69 
Vrsnak, B., Veronig, A. M., Thalmann, J. K., & Zic, T. 2007, 

A&A, 471, 295 

Xu, Z., Lagg, A., Solanki, S., & Liu, Y. 2012, ApJ, 749, 138 
Zhang, Q., Chen, P., Xia, C, & Keppens, R. 2012 



